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1. Introduction 

Taking the AdS/CFT correspondence [1] as a guideline, geometries with an anisotropic 
scaling have been presented in [2] as candidates for gravitational duals for quantum 
critical condensed matter systems that are invariant under Lifshitz scaling 

t -)> A^t, x -)■ Ax , (1.1) 

with a dynamical exponent z > \. A detailed understanding of quantum critical 
metals poses an challenge in theoretical physics [3]. Using gravitational duals to shed 
new light on these and other condensed matter problems continues to be an active 
field of research (for reviews see e.g. [4, 5, 6]). While the validity of this formalism 
is still a matter of debate, it is important to develop these holographic duals further 
in order to be able to test them against experimental results. This paper provides a 
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prescription for defining thermodynamic quantities for dyonic Lifshitz black branes 
that satisfy the expected thermodynamic relations. Such a prescription allows to 
investigate duals of systems with anisotropic scaling in the presence a magnetic field 
- along the same lines as e.g. AdS/CFT duality has been used to obtain a holographic 
description of the Hall conductivity [7]. 

The paper is structured as follows. In section 2, the effective action of the holo- 
graphic model is introduced. It consists of an Einstein-Hilbert term coupled to a 
Proca field and a U{1) gauge field, which is known to give rise to charged Lifshitz 
black brane solutions. Building on previous analysis [8, 9, 10, 11], the model with 
both electric and magnetic fields present is investigated with an emphasis on holo- 
graphic renormalization and thermodynamics. The focus will be on 3+1 dimensions, 
i.e. a strongly interacting field theory dual in 2+1 dimensions. This setup is relevant 
for a holographic description of materials where the charge carriers are confined to 
layers orthogonal to a magnetic field. Another benefit of working in this dimension 
is the duality between magnetic and electric fields, which simplifies the following 
analysis. 

In section 3, counterterms are introduced such that the action and its functional 
derivatives are well-defined and finite on-shell. This extends the analysis in [11], 
where no gauge-field was considered, and allows to renormalize the action for values 
z > 2, which goes beyond the parameter range considered in [12]. The prescription 
presented here is based on a new approach to identify the degrees of freedom of 
the system. Within the space of solutions, Lifshitz-spacetimes form a subset which 
is disconnected from other classes of solutions. Thus, on-shell variations must be 
constrained such that they do not lead away from that subspace. 

The renormalization procedure is then used in section 4 to define the internal 
energy and Helmholtz free energy of the black brane solutions. These obey the same 
relations as in standard thermodynamics. This extends the work of [9, 12] and gives 
further justification that the thermodynamic description of black holes, which is 
known to be valid in the AdS-case, is also applicable for non-relativistic holography. 



2. The holographic model 

The effective action used is of the form introduced in [8] with a Maxwell term added, 

i^bulk ^ 1^ Einstein + JProca + ^Maxwell ■ \^-^) 
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In the above, 



SEinstein — - — T^J / (-R — 2A)wm , (2.2) 

2kL^ Jm 
Sproca = -T^ / (dP A *dP + cP A P) , (2.3) 



S Maxwell = —-. 7^ / FAF, (2.4) 

4/tL^ Jm 

where P is a 1-form on the manifold M and F = dA the Maxwell-tensor. Variation 
with respect to the metric g^y gives the Einstein equations 

G,. + ^g,. = t;, + T^^/ . (2.5) 

The energy tensors T^ and T£/^ are defined as 

T^. = I (p.p. + [dPUdP]/ - \P^P' - \[dPU[dPY-^ , (2.6) 

2 l^^^^^'^ 4 
Variation with respect to P and A leads to the equations of motion 



r.r = \ (f,xF.' - Wf'^ ) . (2.7) 



d*dP = -c*P , (2.8) 

rf*P = 0. (2.9) 

Equations (2.5), (2.8) and (2.9) are known to have asymptotic Lifshitz solutions with 
dynamical exponent z if 

2L2 ' ^ ^ 

c = § , (2.11) 

where d is related to the dimension oi M hj d + 2 = dim M. 

For concreteness, the case dimM = 4, i.e. li = 2, is used in the following. 
This is partly motivated by previous investigations [10, 13, 14], where results so far 
indicated that the qualitative behavior of the system is mainly characterized by the 
ratio d/z rather than d itself. Thus, using general d would merely clutter notation 
without being much more instructive. Beyond that, the duality between electric and 
magnetic field strength in 3 + 1 dimensions allows for some further simplification. 
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2.1 Dyonic black branes 



For the black branes, an ansatz of a static and stationary metric is chosen via the 
tetrad 



eo 
ei 



L- dx , 

r 

L-dy , 
r 

1 

— L — dr 

rg 



(2.12) 
(2.13) 
(2.14) 
(2.15) 



The metric is then given by Qf^^ = rj^^eA^GBu- The tetrad has been introduced for 
later convenience. Furthermore, the orientation on the manifold M is chosen to be 



vm = es A ei A 62 A eo . 
The Proca-field and gauge potential are parametrized as 

Lr'^'^afg^ dt , 



(2.16) 



(2.17) 

A = Lr~'(j)fg'^ dt + LBqx dy . (2.18) 

Here, Bq is the constant field strength of a magnetic field perpendicular to the xy- 
plane. For later convenience, dP and F are parametrized as 

dP = -V2^Lr'''^f dr Adt , (2.19) 

F = Lr-'+^fpo dr Adt + LBq dx A dy , (2.20) 

where the constant po describes the charge density in the system. The above defined 
ansatz will solve the equations of motion (2.5), (2.8) and (2.9) provided the following 
system of first order ODEs holds 



rf = f {z 
9 



1 



% 



1 



rg 

ra 
rb' 

r(f)' 



2(3 + a^) + -IA 



pI + b, 



2 
0^4 



r" + -b' 
2 



-2a 



9' 



zb + a\ A 



pI + b, 



2 
0»,4 



r' + -b' 
2 



26 -2a 



-20 + 



1 
7 



py - A + 



pI + b, 



0^4 



r" + -b' 
2 



(2.21) 
(2.22) 

(2.23) 
(2.24) 
(2.25) 



A straightforward calculation shows that the quantity 

/ / . pI + B?^ . 1 



D. 



,2+2 



-A 



Po + ^0 ^4 



4 



I 2 2 

+a g 



2abg^ 



z¥ - 3^^ 

„2 _|_ r2 

Po 



(2.26) 
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is a first integral. This quantity will prove useful in deriving an equation of state in 
section 4. 

The system (2.21)-(2.25) has the asymptotic Lifshitz fixed point / = foD,g = 
l,a = b = ^ z — 1, = 0. It is however not the only fixed point, as there would also 

be the possibility / ~ r^^^,g = \/ "^^^^^ , a = 6 = = 0, which corresponds to an 
asymptotic AdS solution. The system can be solved in the asymptotic region r — )■ 
by first linearizing around the Lifshitz fixed point, which gives the asymptotic modes 
of the solution, and then iteratively calculating the descendants of these modes. The 
expansion up to the orders that are relevant for this paper can be found in appendix 
A. A crucial result of this expansion (up to a choice of sign) is 



P= IV^^^ + ^l^o, (2.27) 

where ^ is a scalar function that vanishes asymptotically for asymptotic Lifshitz 
solutions. This will play an important role in the calculation of the on-shell variation 
in the next section. 

3. Renormalization 

Using the equations of motion (2.21)-(2.25), the bulk action can be shown to reduce 
to a surface term on-shell 

In the above expression, vqm is the induced volume form on the surface, i.e. vqm = 
ei A 62 A Cq. Plugging in the expansion from appendix A shows that the integrand 
diverges at the boundary as r — )■ 0. Thus, the action needs to be renormalized before 
the standard thermodynamic gauge-gravity dictionary can be applied. This can be 
done as for asymptotic AdS solutions [15], that is by confining the integration domain 
to r > £ and defining a series of counterterms on the boundary dM^, such that the 
limit £ — )■ is well-defined on shell. It should be noted that the following analysis 
is also made under the condition Si = 0, where Si is a coefficient in the expansion 
from appendix A. For z > 2 {z > d) this would anyway be a necessary condition, 
as otherwise there would be a non-renormalizable mode in the solution, while for 
z < 2 {z < d), this choice greatly simplifies the definition of an energy, as it vastly 
reduces the number of required counterterms to cancel all divergent contributions^. 
Furthermore, it will also be assumed that z < 6 {z < d + A). With this, it will be 
sufficient to consider counterterms which are at most quadratic in F and A, which 
simplifies the analysis. At the same time, most known quantum critical systems in 



^Similar conclusions were drawn in [11]. 
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experimental physics have a dynamical exponent in the range 1 < z < 3, which is 
well within the range considered. It can be expected, however, that the final result 
for the thermodynamic quantities in section 4 remains the same for higher values of 
z, but an explicit expression for the renormalized action would involve terms which 
are quartic^ and higher order in F and A. Finally, the normalization /oo = 1 is used 
throughout. This is not a real restriction, as another value would simply correspond 
to a different choice of the time scale, which could be easily be reintroduced by 
multiplying expressions with the appropriate power of /oo. 
The final form of the action is thus given by 

S = Shulk + Sreg + Set ■ (3-2) 

Sreg IS a regulating term necessary to have a well defined variation in the presence 
of a boundary, 

2^1" Sreg = 2 / Kvqm.^ ' f A A *F - — ^ / Q^-A * F . (3.3) 

JdM^ JdMe v^ "~ ^> JdM 

The first term is the usual Gibbons-Hawking term which ensures that the variation 
with respect to the metric is well defined. K is the trace of the intrinsic curvature 
Kfj,u = V(^n;^), which in the given setup can be calculated as K^i, = 2L'§r^t^^-> where 
h^v = g^v — n^riy is the induced metric on dM^. The second term accounts for 
working in a background with fixed charged density, i.e. tegSF, on the boundary 
instead of a fixed chemical potential.^ This choice comes from the consideration that 
the holographic field theoretic problem in mind for this setup is a sample with 2- 
dimensional (semi-) conducting layers in a constant magnetic field with a fixed number 
of dopants rather than a given chemical potential. In the third term, 2^3 is the Lie 
derivative with respect to 63, which can be interpreted as the normal derivative on 
the boundary. This last term is required to cancel a divergent contribution of the 
preceding term that occurs for z > 2. Introducing -k to denote the induced Hodge- 
Star on the boundary and decomposing the Maxwell tensor as 

F' = ie,F , (3.4) 



the regulating term (3.3) can also be written as 



dM ' 



(3.5) 



2KL^Sreg = 2 f Kvqm, - f A A ^F' - — ^ / H^.A A ^F' . (3.6) 

This makes it more manifest that the last two terms correspond to a Legendre trans- 
formation. In the standard dictionary, a bulk gauge field A will source a current J^ 



^There won't be any terms with odd powers as the equations of motion are invariant under 
charge conjugation A -^ —A. 

^This point is explained in e.g. [5]. 
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in the dual field theory on the boundary. In that theory, F^ has an interpretation as 
the response, i.e. {J) ~ |f ~ F^- For z > 2, it can be read off from (A. 5) that the 
mode with Qv will grow faster than the mode involving Q. Thus, the latter needs 
to be canceled and it is actually A + ,^^^, £,e.^A that becomes the source for the dual 
current.^ The Legendre transformation then interchanges the role of A+ ,^^. £,e^A 
as source and F^ as response. 

Finally, Set is a counterterm that cancels all remaining divergences from the bulk 
action. It is given by 

2 2{z + l) f , V2z{z-1) f 

2kL Set = F / ^dM + r / ^^dM 



L JdM L JdM^ 



2(2 -z\ 



[ £e3^ A ^£^3^ + W77^, / F- A ^F™ , (3.7) 

JdM ^[^ — Z) JdM 



where ^ was introduced in (2.27). The first term in (3.7) is simply a boundary cos- 
mo logical constant. The second term cancels a divergence with exponent Z2 coming 
from S2 having a non- vanishing value (cf. appendix A).^ For 1 < z < 2 these two 
terms would actually suffice, for z > 2 there are however further divergences occur- 
ring due to the electric and magnetic fields not falling off fast enough. This is cured 
by adding the terms in the second line of (3.7) involving £e3^ and F"*. 

Two issues might need some clarification. First, the attentive reader may have 
noticed that there is a certain redundancy in the notation as 2esA = F^. This is 
done on purpose to make the conceptual difference between those two terms manifest. 
S^e^A is defined through A, which is the field that enters the bulk action (2.1). 
The field F^ is introduced by performing a Legendre transformation. In the logical 
order, this transformation is done after the action was renormalized. Therefore, the 
counterterm (3.7) is written in terms of Ze^^A and without any explicit dependence 
on F^. 

Second, the first two terms in (3.7) which cancel all divergent contributions when 
A = are different from the terms proposed in [11], where 

2KL^Sa =y f V9M+ ^^""^f" ~ ^^ / V{P,P)VdM . (3.8) 

-^ JdM -^ JdM^ 

However, by using (2.27), a short calculation reveals 

A+^2z{z-l)y^{PJ^) = {z + 3)- "-{P,P) + 0{e) 

= 2{z + 1) + ^/2z{z-l)^ + 0{e) ■ (3.9) 



■^This works analogous to the discussion in [16]. 

^It is worth noting that the ^-term also cancels the divergence proportional to zi when Si 7^ 0. 
To renormalize the action, however, terms with higher powers in ^ would also need to be included. 
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Thus, when ^"^vqm will vanish for r — ;■ 0, which is indeed the case in the parameter 
range considered here, the first line of (3.7) would give exactly the same contribution 
as (3.8). 

Energy and momentum can also be calculated along the lines of [15]. The proce- 
dure differs, however, in the following two ways. First, as was already pointed out in 
[11], the dual theory is not relativistic and thus it is less convenient to work with the 
metric and the stress energy tensor T^'^ = -t^ttt^^ but more useful to work with a 
tetrad and ta, where 

*rA = r]ABT — • (3.10) 

OCb 

Energy, momentum, energy flux and stress are then encoded in the components of 
Ta- a second, and more subtle difference is the way the variations are calculated. 
The aim is to consider gravity duals to systems with anisotropic scaling, but, as was 
noted earlier, in addition to the Lifshitz flxed point there also exists an AdS flxed 
point. In fact, the asymptotic Lifshitz spacetimes form an isolated subset of the 
space of solutions that is disconnected from the subset of asymptotic AdS spacetimes. 
Thus, in the same fashion as the covariant derivative on a surface embedded in M^ 
is basically the derivative of the embedding space constrained to be evaluated on 
curves that do not lead away from the surface, the variations must be constrained to 
'curves' that stay inside the subspace of asymptotic Lifshitz solutions. These 'curves' 
are defined by (2.27), which is a direct consequence of making an ansatz that has 
asymptotic anisotropic scaling. Hence, the variation must be performed under the 
constraint that it is not P, but the scalar ^ which is a degree of freedom of the 
system.^ Using this relation, (3.10) becomes 




2kL^to = 2K -eo- 2Keo + \ — + ^ le^dP 



i{A,F^)eo-{A,eo)F^-{F^,e,)A) 
^ {{2e,A, F')eo - (£e3^, eo)F' - (F^ eo)£e3^) 



2-z 



, 2(^ + 1) , ^2z{z-l) ^ 
+ — ^eo + eeo 



L 

2(2^ z) 



^^r. :T ((^es^, iie3^)eo - 2(£e3^, (^^)^eA) 



+ 2(2^ ((^™, ^")eo + 2F- ■ F^ ■ eo) . (3.11) 



^Making this statement simply based on the expansion (A.1)-(A.5) might appear ad hoc. For the 
purpose of this paper it could just be thought of as a mere working assumption, but an investigation 
of the PDE-system (2.5), (2.8) and (2.9) via a Fefferman- Graham like expansion (cf. [17]) reveals 
that the actual degrees of freedom are not the components of P but are defined via projections to 
60,61,62. These more general results will be reported on elsewhere [18]. 



In the expression above, M ■ u denotes the contraction M^^u^ for a 2-tensor M and a 
l-form u. Ti and T2 are given by similar expressions, but without the terms involving 
^. However, only tq will be relevant when defining the energy in the next section. 

For when considering differentials of thermodynamic quantities later on, it is also 
useful to note the relation 



5F^ 



A + ^£,..4 



(3.12) 



3.1 A note about z 



For -2 = 2, and more generally z = d, the asymptotic expansions (see appendix A) 
become anomalous and contain logarithmic terms. Furthermore (3.6), (3.7) (3.11), 
and (3.12) are not well-defined due to the factor of (2 — z) in the denominator. In 
this special case these need to be modified, 

2KL^Sreg = 2 / Kvqm, - A A i.F^ - L Inr £e3^ A *F% (3.13) 

JdMr JdMr JdM 



, 2iz + l) f , V2z{z-l) f 

2kL Set = f / VdM H f / ^VOM 

^ JdM ^ JdMe 

L f L 

/ lnr£e3^ A7^£e3A + - 

'dM 2 JdM 



+ - lnr£e3^A^£e3^ + 7r / InrF^A^F™, (3.14) 

2 Jfi]\/r 2 



2kL\ = 2K -eo- 2Keo - [ J'^ — ^ + ^ ) i^^dP 

-i{A,F')eo- {A,eo)F' - {F\eo)A) 

-L\nr ((£e3A,F^)eo - (£^3^,60)^^^ - (F^eo)£e3^) 



+ - Inr ((£e3^, £e3^)eo - 2(£e3A, eo)£e3^) 

+ - Inr ((F", F™)eo + 2F™ ■ F" ■ cq) , (3.15) 

2fi:L2-— = -^[A + Llnr£e3A] . (3.16) 

4. Thermodynamics 

A temperature is introduced in the dual theory by considering black brane solutions 
with an event horizon at some finite value r = Tq. At the horizon, f ^ fo,a ^ 
ao, 6 — )■ 60) — ^ 00 while g^ -^ g^l — r/ro). As the form of the equations is invariant 
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under the rescaling r — > \r, Bq — )■ A~^i?o, Po ~^ ^~^P0) the horizon can be assumed to 
be at To = 1. The dependence of the solutions on tq can then be introduced by using 
the rescahng backwards. With this simphfication, the relation between the constants 
at the horizon is given by 

9l = -A - Ibl - p1±^ , (4.1) 

ao = ^ , (4.2) 

9o 

00 = -^ . (4.3) 

9o 

Thermodynamic quantities can now be assigned using the same prescription as for 
the AdS case (see e.g. [4]). The value of k, is associated with the number of flavors A^ 

in the dual Beld theory through 1 = ^. The chemical poteutial ^ maguetic field 

strength b and the charge density q can be read off from the asymptotic expansion 
in appendix A, 

Qu ^ B Q V2NIq .^ ^, 

As -Bo and po enter in a symmetric fashion in the equations of motion (2.21)-(2.25), 
the magnetization density is given by 

m= = --^ . (4.5) 

2kL 4/t2 q ^ ^ 

This relation will become more clear in the discussion in subsection 4.2. Reintroduc- 
ing the scahng in vq reveals that m, /x oc r^^ and the values of b and q are related to 
the variables at the horizon via 

h-^ q= ^" (4 6) 



L^r-Q 2KL^rQ 



A temperature is defined via Wick rotating time and then compactifying on the 
thermal circle, the result is 

rr fo9o 



The value of tq also defines the entropy density, 

27r 



kL'^7 



(4.7) 



(4i 



The evaluation of the conserved quantity (2.26) at the horizon and in the region 
r — )■ relates T and s with the variables in appendix A. 

2kL^sT = i ^ l« + t^ j Q 2; ^ z , . . 



10- 



4.1 Thermodynamic potential and equation of state 

By the standard prescription, the grand canonical potential is associated with the 
the value of the renormalized Euclidean on-shell action. However, in the case at hand 
(3.2) contains the term — J^^ I A + ,^^, £,e.jA j A •kF'^ which has been added to the 
action to allow for a setup with fixed charge density. This terms is not part of the 
renormalization to cancel divergences, but it changes the thermodynamic potential 
by the value /iq, resulting in the canonical ensemble. Thus, from the on shell value 
of the action (3.2), 

•\i rji qEucl, on— shell /a M^\ 

V is the volume of the system and a is the Helmholtz free energy density. Plugging 
in the parametrization presented in section 2.1 and using the asymptotic expansion 
from appendix A leads to 

An internal energy can be defined by working in the spirit of the AdS/CFT corre- 
spondence and considering the on-shell action (3.2) as a generating functional for 
the dual field theory with the boundary values of the fields interpreted as sources 
for their dual operators. In [15] the dual stress energy tensor T^"^ was considered as 
the operator that is sourced by the boundary metric /i^i^. As already indicated in 
section 3, instead of T^^", the quantities ta defined in (3.10) will now be used for this 
purpose. With dt being the Killing vector that generates time translation invariance 
the internal energy density e is associated to tq, given in (3.11), through 

Lt = ro{dt) . (4.12) 

sources=0 

The subscript sources = reminds of the fact that according to the standard de- 
scription, the right hand side of (4.12), which comes from a functional derivative of 
the on-shell action, must be evaluated with all sources, i.e. independent boundary 
values, set equal to zero. To account for this, To{dt) must be evaluated at the point 
where the explicit dependence on the source F^ is set equal to zero. Formulated 
quantitatively, this means that the term — J^^ (A + ,^^. £,^3^) A •kF'^ in (3.6) will 
not contribute to the internal energy. This is also sensible, as this term would give a 
contribution coming from having a nonzero chemical potential in the system, whereas 
the internal energy by definition should just account for the mass of the black brane 
that causes the curvature of spacetime. The result is 

4{2-2)Vi^A4 ^ / 2 



2KL-'t = { ' ^ ' (4.13) 

^ -2M z = 2. ^ ^ 
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Now (4.9), (4.11) and (4.13) can be combined to 



a=c-sT 



(4.14) 



which is indeed the correct expression for the density of the Helmholtz free energy. 
Furthermore, from the above calculations an equation of state can be derived. 



z + 2 



sT — mb + fiq 
sT ~rab + nq + 



b'^ + fq'^ 



Zy^2, 

z = 2 . 



(4.15) 



The first line is in accord with the findings in [19] and [12]. The appearance of b^ 
and q^ in the equation of state for 2; = 2 is an artifact of an ambiguity in defining 
a counterterm for this particular value of the dynamical critical exponent. The 
approach presented in section 3 was a minimal one, i.e. just taking the counterterms 
which are required to cancel all divergences. This results in the coefficients in (4.15). 
As a matter of fact, for z = 2 it would be possible to add the terms 2 J^^ -Cea^ A 
T*rF^ — Jg^j £e3^ A -ki^e^A and Jqj^ F™- A -kF"^ with arbitrary coefficients to the action. 
This would leave (4.14) unchanged, but would alter the coefficients of b^ and q^ in 
(4.15). In particular, it would be possible to cancel these coefficients, making the 
second line of (4.15) identical to the first. It is unclear, what argument should be 
used to single out this choice and fix the ambiguity. 

4.2 The differential of the Helmholtz free energy 

In thermodynamics, a satisfies 



da 



-sdT — mdb + fxdq . 



This corresponds to the three relations 

da 
df 

da 
db 
da 
dq 



b.q 



T,q 



-m 



fi 



(4.16) 

(4.17) 
(4.18) 
(4.19) 



T,b 



These can easily be verified for dyonic AdS black branes in the case oi z = 1, where 
an exact solution is known. What will be shown in the following is that they also 
hold for z > 1. 

First of all, (4.19) is a direct consequence of (3.12) when taking the limit r — )■ 0. 
From this, (4.16) will follow if it can be shown that any of the relations (4.17)-(4.19) 
implies the other two. To proceed with the proof of this, it is useful to note that in 
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the equations of motion (2.21)-(2.25) as well as (4.11), the values of Bq and po only 
occur in the combination 

V = B'o+pI. (4.20) 

Furthermore, as the dependence of tq just enters in the form of a rescaling of the final 
expression, all so far introduced thermodynamic quantities must be of the form — ^ 
with s being some scaling exponent and Q a function of a single variable.^ Therefore, 
let the functions J-' and Q be defined via 

T = ^^ , (4.21) 



2k q b r-Q 



2«T = ^^ . (4.22) 



Imposing the conditions ro = 1 and Si = on the ODE system (2.21)-(2.25) results 
in a one-parameter family of solutions, the parameter being 77. Hence, J-' and Q are 
not independent and must satisfy a non-trivial relation. This relation turns out to 
be 

AT' -{z- 2)g + Ar]g' = . (4.23) 

The validity of this will follow as a corollary to what will be proved in the following, 
namely that (4.23) is equivalent to each of (4.17)-(4.19). 
First of all, the differentials of T, b and q are 

zj^ 2T' 

dT = - dr-o + -r-r^iBo dB^ + po dpo) , (4.24) 

^^ = -n^dro + j^dBo , (4.25) 

From this, at constant b and q, 

1 



dT 



{-zJ^ + Avr)dro, (4.27) 




^(5T)[^ = ^^^^[-(z + 2)^ + 4r^^V^o, (4.28 



d{^ 

q 



9k-/' 

= ^n[-(^-2)^ + 4r7^V^o. (4.29) 

b,q ^0 



^Of course, q, /i, b and m are not exactly of this form, they however differ only by a factor of q 
or b respectively. 
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Therefore, 
dT 

fa,q 



d_ 
df 



T^S -h 

2 + 2 z + 2 4/t2 q 



b,q 



27r 



KLVg {zj= - Ar]r) 



zj= 



4:Z _, 2(^-2) ^ 



2r" 



z + 2 



z + 2 



z + 2 



v'Q^ 



27r 47rr7[4J^' -{z-2)g + 4r]g'] 



kL'^tI {z + 2)KL'^rl{zJ^ -AriP) 

—5 . 



(4.30) 



The last equahty follows from (4.23). This establishes the equivalence of (4.23) and 
(4.17). 

In an analogous way for constant T and q, 

1 



d 



db 

ds 



T,q 2L^Borlr 

4:71 , 



- {zj^ — 4:1] J^') dr-Q 



Ta 



Ta 



K.L'^r, 



(4.31) 
(4.32) 







2kI 

[-(z-2)rg + zTg']dro 



r-o '^' 



kL 



2rl-'r]J^' 
In the last line, (4.23) was inserted. With this. 



- [{z - 2\zT - 4r]J^')g - 4zJ^J^'] dr^ 



(4.33) 



da. 



ta 



d_ 



z ^ 2 b2 + 4/t2q2^"i 

z + 2 z + 2 4/t2 q 



z ^ds 



2 



1 bfi 2 b2 + 4/t2q2^ d f fi 



z + 2 db ^^^ ' z + 2 n'^q ' z + 2 
1 b/i z-2LBo{b^ + 4K^q^ 



4k^ q (9b \ q 



ta 



+ 



z + 2H?q z + 2 
Borl 



g 



+ 

b^ 
4fi:2q 
— m . 



Lk {zJ= - 4^J=') 



4/tro '^r] 
■ 2z Tr 2z L4(b2 + 4/t2q2)jrjr' 



z + 2 rg+3 z + 2 



rl ^r] 



(4.34) 



Again, as the equality holds if and only if (4.23) is assumed, the equivalence of that 
assumption to (4.18) is proved. Due to the symmetric appearance of b and q, this 
must also be true for (4.19). As the validity of (4.19) has already been established, 
this concludes the proof of (4.16). 

Unfortunately, the thermodynamic relations found so far are not sufficient to 
determine an explicit expression for J-" or ^. A few exact solutions are known (see 
appendix B), but in general numerical methods are needed to study these functions. 
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4.3 Susceptibility and magnetization 

It is also possible to derive an expression for the density of the magnetic susceptibility, 
dm 



X 



db 



m 



'^bQ 



T,c\ 






Anr'r^ 



V 



{z - 2)g - 



2(b2 + 4/t2q2) 



2„2 



[zb^ + Sft^q 



zF -ArfF' _ 
m zLb'^T 



z'Lb'r 



2 772 



AkTq "^1] 



Anr^ ^r]{zT - 4r]T') 



Using the specific heat at constant volume, 



Cv 

T 



ds 
df 



IQtt^ 



kLtI-' 



as well as (4.8) and (4.21), (4.35) could also be written as 

[zb"^ + 8fi:2q2) tn zTb'^{2s - zcvT) 



X 



2(b2 + 4/t2q2) b 4(b2 + 4/t2q2) 



(4.35) 



(4.36) 



(4.37) 



As a consequence, in the limit of vanishing temperature, assuming s and c\;T do not 

diverge in this limit*, 

(zb2 + 8fi:2q2) m 



^T=o 2(b2 + 4fi;2q2) b 
and for vanishing magnetic field, 

m 1 /i 



X 



lim ^ 

b-^O b 



(4.38) 



(4.39) 



b=o b-^o D 4/t2 q 

For z = 1 this identity can easily be checked for dyonic AdS black branes. That it 
also holds for 2; > 1 based on (4.16) is a non-trivial result. 

As was noted in [12], for 1 < z < 2, the value of z/ from appendix A can be 
expressed as 



u 



^z-l 



fdr . 



(4.40) 



Thus, because / > outside the horizon, from (4.5) follows that m has the opposite 
sign to b. As this also implies that x will be negative, at least in the limits of 
low temperature and magnetic field strength, this means that the system exhibits 
diamagnetic behavior. In contrast, for 2 < 2; < 6, an expression for u is given by 



u 



r'-\f-l)dr- 







1 

2-2 



z = 2, 
2< z <6 . 



(4.41) 



^Numerical investigations in [13] suggest that they remain finite. 
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Here it is potentially possible to have a setup with m and b having the same sign 
and X positive and thus modeling a paramagnetic material. In fact, the known exact 
solution for 2; = 4 (see appendix B.2) is such a case. 

4.4 Numerical results 

Though the main results of this paper are derived analytically, it is instructive to also 
have a numerical check of certain equalities. As the qualitative features seemed to 
be rather indifferent to the particular value of z, numerical results are just presented 
for one value, 2: = |. 

The information about thermodynamic quantities is encoded in the functions J-" 
and Q (4.21) and (4.22). A plot of these functions can be seen in figure 1. The 




1.95 




Figure 1: Plot of the functions J-{i]) (left) and Q{ri) (right). F vanishes for r] = ^^, 
indicating that the black brane becomes extremal for this value. 



function J-" can be seen to vanish sublinearly at the critical value 77 = — 4A = ^. 
As this happens at a finite value of tq, this means that the black brane becomes 
extremal. A more detailed numerical investigation of the behavior of the system 
when approaching criticality can be found in [10]. 

The aim now is to check (4.23). To do this, define the functional 



J^ = J^G-r]G + 



4 



(4.42) 



which is the general solution of (4.23) for given Q, and fix J-'o such that J-" and J-' 
coincide at some value of rj. Then (4.42) can be compared with the numerical value 
of J-" at other values of rj. The relative error 



^rel^ 



7 -T 



7^7 



(4.43) 



is plotted in figure 2. For a better comparison later plots, t] has been translated 
back into a value of temperature, normalized by the temperature at q = 2^^- For 
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0.001 



0.01 



Figure 2: The relative error AreiJ^ versus temperature. The loss of precision for lower 
temperatures is due to the numerical function values become less precise in that region and 
due to an accumulation of numerical errors from the integration. 



temperatures of 0(1) and higher, the deviation can be seen to be lower than ten 
significant digits. When the temperature is lowered, the deviation increases. This can 
be attributed on the one hand to numerical values of J-" and Q having lower precision 
at low temperatures and on the other hand to the accumulation of numerical errors 
when integrating (4.42). 

As (4.23) was shown to be equivalent to (4.16), the above results give a good 
indication that the latter is indeed satisfied. It is however also be possible to make 
a more direct check. For this purpose, the relative error 



Aj-elS 



da\ 
dTlb, 



+ 3 



da\ 
9Tlfa,q 



— 5 



(4.44) 



is plotted in figure 3. Results here are less precise than for (4.43). This is mainly due 
to the fact that a numerical estimate of a derivative is in general more susceptible 
to the precision of the input data than the estimate of an integral. Nevertheless, 
the numerics show an agreement to at least three significant digits - and even better 
agreement for higher temperatures, where the precision is better. 

This section concludes with a numerical check of (4.39), which also was a con- 
sequence of (4.23). This identity allows to compare a second derivative of a with 
quantities that can be read off from the asymptotics. Figure 4 shows the relative 



error 



Are/X 






T,q 



1 M 

4k2 q 



ab^lT.q 



J 1_M 

"^ 4k2 q 



(4.45) 
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0.001 



10- 



10- 



10- 




10- 



0.001 



0.01 



0.1 



Figure 3: The relative error As versus temperature. The increase for fower values of 
temperature is due to the lower precision of the numerical data in that region. 



at vanishing magnetic field. Also here an agreement of about three significant digits 



l^retX 



0.001 



10- 



10-7 



10- 



..^^ 






10- 



0.001 



0.01 



0.1 



Figure 4: The relative error AreiX versus temperature for b = 0, normalized with the 
values at q = 2^. As in the other plots, the error has a tendency to increases with lower 
values of T due to precision issues in that region. 



or better can be seen. 

The results presented here were as far as it was possible to go using numerics 
with reasonable computation time. The general trend was that the deviations pre- 
sented above decreased when the precision was increased. It stands to reason that 
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computations done with an even higher number of significant digits would further 
improve the numerical results. 



5. Conclusions 

In this paper, several results on of dyonic Lifshitz black branes were established. 

Though the renormalization of the action presented in this paper is still work in 
progress, preliminary results for black branes with d = 2 and z < 6 can be obtained. 
The task of deriving an expression for general values of d, higher values of z and even 
a curved horizon seems straightforward, however, an implementation is expected to 
become more complicated when expansions up to higher order need to be taken into 
account. 

A further result would be the evidence that this renormalization procedure gives 
expressions for the Helmholtz free energy and internal energy that are in agreement 
with standard thermodynamic relations. This bolsters the case for dyonic Lifshitz 
black branes as candidates for a holographic description of phenomena involving 
magnetism in quantum critical systems with a non-trivial dynamical exponent. 

Finally, the magnetization and susceptibility of the dual theories were worked 
out using the gravitational description. For 1 < z < 2, the low temperature limit 
is always diamagnetic, whereas paramagnetism can occur for 2 < z < 6. The 
only known exact solution for z = 4 happens to be of the paramagnetic type, but 
it remains an open question whether paramagnetism is the rule in this region of 
parameter space. 
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A. Asymptotic Expansions 

An expansion for generic z for asymptotically Lifshitz solutions of (2.21)-(2.25) is as 
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follows. 



f = fo 



1- 



2^/^^l{z, - 2)E, ,, ^/^^1{Z2 - 2)E2 



^zi _ 



,Z2 



z{-2 + 2z + zi) z{-2 + 2z + Z2) 



2 + z 



z (^2 + 3^ + 2) 



„2+z 



[z - l)iz - A){B' + Q'l^, ^ ^ ^ 



9 = 1 



8(^-2)2(2 + 1) 

y/z - IZiEi 



r"i + 



y/z - \z{i 



+ 



z(-2 + 2z + Zi) ^(-2 + 2^ + ^2) 

(223_5^2^3^_2)SiS2' 



V-2 - 1-M + 



2^2(^ + 1] 



„2+z 



(2z - 3) (g2 + Q2) ^ 
4(z-2)2(z + l) 

zi(zi -2)Si 



^(-2 + 2^ + ^1) z(-2 + 2z + Z2) 

(^-4)v/^^T(i32 + Q2^ 



V-2- 1- 



4(^-2)2(^ + 1) 

2r^i2iSi 



-r^ + 2A^r2+^ + . . . 
2r^2^2H2 



2(-2 + 2z + 2i) ^(-2 + 2^ + 22) 
(^_ 4)7^^1(^2 ^g2 



^'\4 4A^2+. 



4(^-2)2(^ + 1) 



-r — 



-r 



2; 



Q 



-r''z/ + 



-r^ + . . . 



(A.i; 



(A.2) 



(A.3) 



(A.4) 
(A.5) 



where . . . indicate descendants of the previous listed modes. The Exponents z\ and 
Z2 are given by 



Z\ 



Z2 



z + 2-^{2 + zf + 8(^ - l){z - 2) 
z + 2 + V(2 + z)2 + 8(z - l){z - 2) 



(A.6) 

(A.7) 



For z > 2, the exponent Zi > and thus, in order to be renormalizable, the solution 
must have Si = 0. For the marginal case z = 2 with zi = Q and Z2 = 4, there 
are logarithmic modes occurring. After discarding a growing mode, the expansions 
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(A.1)-(A.5) are modified to 



f = l + r' 



32M - 32E2 - B' - Q' 8M + B^ + Q 



+ 



96 



16 



Inr 



16 



(Inr) 



(Ai 



g = l+r^ 



2 I ^2^, 



28M + 32E2 + B' + Q' -l6M + 3iB' + Q 



16 



96 
-(Inr) 



32 



Inr 



B' + Q\^ ,2 



+ ... , 



(A.9) 



a = 1 + r^ 



2 I n2\ 



M-6AE2 + B' + Q' 1QM + 3{B' + Q 



96 



16 



Inr 



2 I ^2 



B' + Q 



(Inr) 



6 = 1 + r^ 



20M - 32E2 - B' - Q' 16M + B' + Q 



2 I n2 



+ 



48 

B^ + Q^ 



16 



Inr 



(Inr) 



+ ... 



Q [r^(-z/ + lnr) + . . .] . 



(A.IO) 



(A.ll) 
(A.12) 



B. Exact solutions 



B.l z=l : dyonic AdS black holes 



The solution of a dyonic black brane with a horizon at r = 1 is given by 



/ = l,/ = l-(l + ^^)r^ + ^^r^a = 6 = O,0=p(-r + r^) . (B.l) 
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From this follow the thermodynamic quantities 



Po 



2/tL2r2 ' 






2n 



kL^Tq 



2KLh^ 
-4 + 3p2 + Wl 



cv 647r2 



) 



T SkLto (4 + p^ + ^2 

.r-o (12 + 3pg + i 
6«:(4 + p2 + 52 



2\ ' 



B.2 z=4 



For 2; = 4, there is the special solution^ 



(B.2 



p = - — , [D.S 



(B.4 



— -^0 m p: 

2KLro 
T = '-^ ^ , B.6 

1671 Ltq 



(B.7 



e=i±4+Jl. (B.8 



(B.9 
(B.IO 



Lro (12 + 3pg + i?g) 



/ = 1,/ = 1 - r^a = v^,6 = v^,<l> = - ^/f' . , (B.12) 

2 (r^ + 1) 



which solves (2.21)-(2.25) provided r] = B^+pl = 8. As it is just an isolated solution 
at a single value of rj, the thermodynamic quantities that involve differentiation can 



^This basically is the solution presented in [10], rotated on the -BqPo plane. 
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be calculated. The ones obtainable 


are 


q = 


Po 


2/tL2r2 ' 




Po 


/i = 


2Lrl ' 




^0 


b = 


LVg ' 


m = 


Bo 


AnL'^rl ' 




1 


T = 




TrLr^ ' 


5 = 


27r 


/tL^rQ 


e = 


0, 


a = 


-2 



kL^i 



(B.13) 
(B.14) 

(B.15) 

(B.16) 

(B.17) 

(B.18) 
(B.19) 
(B.20) 



It might be worth noting that in this solution q and /x have opposite sign and the 
internal energy is vanishing. 
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